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We present the single stellar population (SSP) synthesis results of our new syn- 
thetic stellar atmosphere models library with a spectral sampling of 0.3 A, covering 
the wavelength range from 3000 A to 7000 A for a wide range of metallicities (twice 
' solar, solar, half solar and 1/10 solar). The stellar library is composed of 1650 spec- 

, tra computed with the latest improvements in stellar atmospheres. In particular, it 

■ incorporates non-LTE line-blanketed models for hot (Tcff ^ 27500 K), and LTE line- 

' blanketed models (Phoenix) for cool (3000 ^ Tcs ^ 4500 K) stars. Because of the high 

l/^ , spectral resolution of this library, evolutionary synthesis models can be used to predict 

' the strength of numerous weak absorption lines and the evolution of the profiles of 

"^H I the strongest lines over a wide range of ages. The SSP results have been calculated for 

Q^i ages 1 Myr to 17 Gyr using the stellar evolutionary tracks provided by the Geneva and 

I ' Padova groups. For young stellar populations, our results have a very detailed coverage 

O I of high-temperature stars with similar results for the Padova and Geneva isochrones. 

For intermediate and old stellar populations, our results, once degraded to a lower 
^ ' resolution, are similar to the ones obtained by other groups (limitations imposed by 

the stellar evolutionary physics notwidthstanding) . The limitations and advantages of 
our models for the analysis of integrated populations are described. The full set of the 
stellar library and the evolutionary models are available for retrieval at the we bsites 
^ . fhttp : //www, iaa. csic ■ es/'-^rosa and http: //www, iaa. csic . es/^mcs/Sed@( or on 

' request from the first two authors. 
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1 INTRODUCTION 

Evolutionary synthesis, introduced by iTinslevI il96d) . is 
a powerful tool for interpreting the integrated light of a 
stellar population. This technique makes a prediction for 
the properties of a stellar population taking as a free pa- 
rameter the star formation rate, the initial mass func- 
tion (IMF) and the metallicity. During the past decade, 
evolutionary synthesis models have been developed and 
improved significantly, using the newest and more ad- 
vanced stellar evolutionary tracks and stellar libraries 
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they have been exte n ded to new wavelength domains (e.g . 
ICervino et al] l2000l : ICervino. Mas-Hesse, fc KunthI l200l 
for 7 and X-rays respectively). However, it is also well 
known that the power of this technique is very dependent 
on the reliabilit y of the stellar models and libraries used 
as in gredients iLeitherer. Fritze- v.Alvensleben. &: Huchral 
Il996l and references therein), as well as on the 
size of the system to which the models are ap- 
plied l|C ervino. Luridiana. fc Cast andeJ200otlBruz uall2002| : 
Cervino et al.l l200ll 120021: ICeryino fc Luridianal l20o3 : 
Gilfanov. Grimm, fc SunvaevI l2004ir One important ques- 
tion related to this technique is whether it is actually up to 
the challenge of interpreting multi-wavelength data to infer 
the properties of stell ar populations of v e ry different ages 
and metallicities (e.g. lOriglia et al.lll999l : iMassev fc OlsenI 
|2003). 

In recent years, numerous observational surveys have 
been undertaken to measure the redshift evolution of the 
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stellar properties of distant galaxies, with the goal of under- 
standing the formation and evolution of galaxies. Spectral 
observations of intermediate and high redshift galaxies have 
shown a large variety of stellar populations. The integrated 
light of these galaxies has been analyzed with evolution- 
ary synthesis models to derive ages, masses and metallici- 
ties. At ultraviolet (UV) wavelengths, many of these galaxies 
show the footprints of recent bursts of star formation. The y 
are classified as star forming galaxies ijSteidel et alj[l99^ . 
The UV spectra of these galaxies, which are very similar to 
local starbursts, are dominated by resonance lines formed 
in th e interstellar mediu m and/or in the winds of massive 
stars jPettini et al .'2000V Observations at rest- frame optical 
wavelengths of these high redshift galaxies (z^ 2) are scarce 
due to the lack of non-cryogenic near-infrared (NIR) multi- 
object spectrographs in very large telescopes. Ha nebular 
emission has been detected in a few of t hese galaxies. How- 
ever, as in nearby star forming galaxies l)Brinchmann et all 
|2004) . the higher order Balmer lines appear in absorption, or 
in emission with strong absorption wings. Thus, to estimate 
the star formation rate in these galaxies requires a careful fit 
of the stellar continuum and subtraction of the absorption 
component of the Balmer lines. Evolutionary synthesis mod- 
els at intermediate or at high spectral resolution are needed 
to perform this fit. 

With this motivation, Gonzalez Delgado, Leitherer, & 
Heckman (1999, hereafter GLH99) built synthetic spectra of 
H Balmer and He absorption lines of starbursts and post- 
starburst galaxies, usi ng as input to the evolutionary code 
Starburst99 (Lcithcrcr et al.lll999l) a stellar library of syn- 
thetic spectra, with a sampling of 0.3 A, effective tempera- 
ture from 50000 to 4000 K, a nd gravity in the range 0.0 ^log 
g^ 5 .0 at solar metallicity iGonzalez Delgado fc Leitherer! 
I1999I . hereafter GL99). The line profiles were computed us- 
ing a set of programs developed b y Hubcny and coUabora- 
tors llHubenvl[l985 iHubenv. Lanz. Jeffervlll99^ . Due to 
computational restrictions, we only computed small spectral 
ranges around the most important Balmer and Hel lines. 

Many beautiful data sets on starbursts and 
galaxies with active galactic nuclei (AGN) have 
been obtained at intermediate spectral resolution 
in the last few years to date the stellar popula- 
tion 
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These results have motivated us to build a much improved 
stellar library, both in the mapping of the {T^s, log g) 
plane and in the complete coverage of the optical spectral 
range at high spectral resolution. With this new library 
it is possible to predict the spectral evolution not only 
of young starbursts but also the properties of old stellar 
systems. Intermediate or high spectral resolution mod- 
els are also required to constrain the stellar population 
in el lipticals, and to break the age-metallicity degener- 
acv iVazdekis fc ArimoTol Il999f) . With this motivation, 
iBruzual^^Cha^o^ l200a) presented models, computed 
with their code GALAXEV, at intermediate resolution (3 
A) covering a large range in ages (0.1 Myr to 20 Ga), 
and wavelengths (3200-9500 A) for several metallicities. 
GALAXEV models are in some aspects similar to those 
presented in this paper, but they have much lower spectral 



resolution, different coverage of the (Toft , logg) plane and, 
most importantly, they use empirical stellar spectra. A 
comparison of the two sets of models is presented Section 5. 

This paper is organized as follows. In Section 2, we 
briefly discuss the stellar library. A more detailed descrip- 
tion is given bv lMartins et al] ||2004) . Section 3 describes the 
main stellar ingredients, the assumptions, and the computa- 
tional techniques of the synthesis models. Section 4 presents 
the results and a discussion of the spectral evolution of stel- 
lar populations, at several metallicities, which evolve accord- 
ing to the Geneva or Padova isochrones. The models are 
compared with previous work in Section 5, and with obser- 
vations of star clusters in the Large and Small Magellanic 
Clouds in Section 6. We discuss the limitation of the models 
in the red supergiant phase in Section 7. The summary and 
conclusions are in Section 8. 



2 STELLAR LIBRARY 

This stellar library is an extension of the library built by 
GL99. Her e, we present a brief description. A full discussion 
is given in lMartins et alj i2004ft . 

The grid includes the synthetic stellar spectra from 3000 
to 7000 A with a final spectral sampling of 0.3 A. The spectra 
span a range of effective temperature from 3000 to 55000 K, 
with variable steps from 500 to 2500 K, and a surface gravity 
logg = —0.5 to 5.5 with dex steps of 0.25 and 0.5. For each 
temperature, the minimum gravity is set by the Eddington 
limit. The library covers several metallicities: twice solar, so- 
lar, half and 1/10 solar. We assume solar abundance ratios 
for all the elements, and a helium abundance of He/H=0.1 
by number. In addition to the spectral and metallicity cov- 
erage, this new library presents the following improvements 
with respect to GL99: 

(i) Hot stars. Line-blanketed, non-LTE, plane-parallel, 
hydrostatic atmospheres are used to compute the stellar 
spectra of O-type star s. The spectra are from the grid of 
iLanz fc Hubenvl ^|2003^ ^. This grid has 12 effective temper- 
atures, 27500 ^ Teff ^ 55000 K, with 2500 K steps, and eight 
values of the surface gravity, 3.0 ^ logg ^ 4.75, with 0.25 
dex steps. These models assume a microturbulence velocity 
of 10 km s"\ 

(ii) Cool stars. Line-blanketed, LTE, spherical at- 
mospheres are computed with the PH OENIX code 
jHauschildt fc Baronll999l:lAllard et al.l200J) for stars with 
3000 =^ Teff ^ 4500 K, with steps of 500 K. The grid covers 
13 values of surface gravity, —0.5 ^ log <? ^ 5.5, with 0.5 dex 
steps. The models assume a mass of 1 M0, a mixing-length 
equal to 2.0, and a microturbulence velocity of 2 km . 

(iii) Intermediate temperature stars. As in GL99, the syn- 
thetic spectra for stars with 8000 ^ Tcff ^ 27000 K ar e ob- 
tained with Kurucz LTE a tmo spheres'^ (Kurucz||l9 93() an d 
the program SYNSPEC^ llHubenv. Lanz. fc Jeffervl Il995^ ■ 
except that for stars with 4750 ^ TcS 7750 K, we 
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use the program SPECTRUM* iCrav fc Corballvl Il994) . 
together with Kurucz atmospheres. Differences found be- 
tween the metailic line strengths and molecul ar bands using 
SPEC TRUM and SYNSPEC are discussed in lMartins etaUi 
J2004I) . 

In all the cases, the maximum distance between two 
neighboring frequency points for evaluating the spectrum is 
0.01 A. However, the spectra are degraded to have a final 
resolution of 0.3 A by performing a rotational and instru- 
mental convolution for each spectrum. 

The library contains 415, 409, 416 and 411 spectra for 
Z=0.002, 0.010, 0.020 (solar) and 0.040, respectively. The 
typical grid coverage is illustrated in Figure 1. Isochrones 
obtained with the evolutionary tracks from the Geneva and 
Padova groups (cf. Section 3) have been over-plotted. The 
figure illustrates that this library constitutes a homogeneous 
set of stellar spectra covering a dense grid of fundamental 
parameters, including a significant number of models that 
reduce the problems of the close model as signation, assumed 
in most synthesis codes (see, however, I.Tamet et alJl20o3. 
for a different approach). This library al so presents some 
impr ovements with respect to STELIB^ llLe Borgne et all 
12003^ . the stellar library used by GALAXEV; in particular, 
it includes models with T^a > 25000 K, and T^tt 4000 
K and logg 1.0, which allow the prediction of the spec- 
tral evolution of the stellar populations of both very young 
starbursts and old systems. Furthermore, the high spectral 
resolution of the library allows us to predict the evolution of 
the Hei and Hen lines, and to estimate the absorption cor- 
rection to the nebular He emission lines. This correction is 
very important for an accurate determination of primordial 
He abundance using nebular emission lines of star forming 
galaxies. 

The full set of the stellar library is available for retrieval 
at the website ( http: / / www.iaa.csic.es/~rosa|f or on request 
from the authors. 



3 DESCRIPTION OF THE INPUTS 

The library has been implemented in Starburst99 and 
Sed@^. There are no significant differences between the re- 
sults obtained with the two codes when models with the 
same ingredients (evolutionary tracks, IMF, ....) are com- 
pared. Here, we describe only the assumptions and prescrip- 
tions related to Sed@ models. 

There is, however, a principal difference in the way the 
two codes. Starburst99 computes the isochrones from the im- 
plemented evolutionary tracks and obtains the results from 

* Available at 

http : //wwwl . appstate . edu/dept/physics/spectrum/spectrum.html 
" Available at http://webast.ast.obs-mip.fr/stelib/ 
^ Sed@ is a synthesis code included in the Legacy Tool project 
of the Violent Star Formation European Network. The code is 
written in ANSI C under GNU Public License and, currently, is 
coordinated by M. Cervino; the code and their results must be 
referred to solely by its documentation, Sed@ Reference Man- 
ual (in preparation), its web address and the citations in the 
headers of the resulting files. More information can be found at 
[http : //www. iaa. csic . es/~mcs/Sed£J 



the implemented atmosphere library. In contrast, Sed@ cur- 
rently does not compute isochrones, but takes them from 
the literature if they are in the appropriate format. 

3.1 Stellar initial mass function (IMF) and star 
formation law 

The results presented here assume a power law IMF, $(m) — 
Cm", where the constant C is determined by the total gas 
mass co nverted into sta rs. The slope a — —2.35 is that de- 
rived bv lSalpeteil 1^5^. The low-mass and high- mass cut- 
offs have been ch osen as Miow =0.1 M p, and Mup ~ 120 M0 , 
respectively (see lGirardi fc Bicalll99j. for a detailed discus- 
sion about the more appropriate choice of the lower mass 
limit). The IMF is binned into a grid of the initial masses 
determined by the particular isochrone, and all the stars be- 
longing to the same mass bin are assumed to have exactly 
the same properties. The synthesis is done directly by adding 
the different contributions from the stars to all the synthe- 
sised wavelengths. Here, we explore only instantaneous burst 
models, normalized to IMq. 

3.2 Stellar evolution and isochrone integration 

The models have been computed using the isochrone tech- 
nique to produce a sm ooth spectral evolut i on. T his tech- 
nique was introduced bv lCharlot fc Bruzuall il99ll) to avoid 
abrupt changes of the properties of the stellar popula- 
tions associated with changes in the evolutionary stages 
of the stars when the mass resolution of the stellar tracks 
is inadequate. The isochro ne synthesis techn i que is differ- 
ent from the one used by iTinslev fc Gumj (|197d), takes 
into account the lifetime of the stars in each evolutionary 
phases considered, imposing that the number of stars in a 
given phase of the stellar evolution must be proportional 
to the life time of such phase of an individual star. Such a 
condition, the so ca lled Fuel Consumption Theorem, FCT 
iRenzini fc Buzzonilil986i : iBuzzoni«il989i') . is specially rele- 
vant for the faster phases of the stellar evolution. Although 
both methods should be equivalent, the results depend on 
the method used (see Buzzoni 1995, 2002b, for an extensive 
discussion) . 

In this work, two different sets of stellar isochrones are 
used. Th e first one is the set of is ochrones from the Pa dova 
group by'Bertelli et al' ('1994') and lGirardi et al.l ll2000h pre- 
sented in Girardi ot al. (2002)'^, with metallicities Z=0.019 
(solar), 0.008 and 0.004 (hereafter Padova isochrones). The 
initial mass resolution of the used Padova isochrones is suf- 
ficient to assure a correct integration of the isochrone-IMF 
product. These models follow the evolution from the zero- 
age main sequence to the early asymptotic giant branch 
(AGB) and to the first thermal pulse, or to the carbon 
ignition depending on the initial mass. They cover an age 
range from 4 Myr to 17 Ga in a set of 74 isochrones with 
a 0.06 log time sampling. An additional set at Z= 0.030 
and 0.002 is also available on the Padova web server but 
it comprises a smaller range of ages (0.063 to 15.8 Gyr). 
We note that there are different sets of isochrones released 
by the Padova group, in particular, the isochrones from 

^ The full set is available at |http : //pleiadi .pd. astro ■ it/| 



4 R.M. Gonzalez Delgado et al. 




log Trff (K) ^ l„g Teft (K) 

Figure 1. HR diagram: Points represent the stellar library grid. Lines are isochrones of stellar populations that evolve following the 
Geneva (left) and the Padova (right) evolutionary tracks at metallicities Z=0.040 and Z=0.008, respectively. The isochrone ages are from 
1 Myr to 10 Gyr for the Geneva and 4 Myr to 17 Gyr for the Padova models. 



iMarigQ fc Cxirardil (|200J^ fulfill the FCT requirements, un- 
fortunately sucli isochrones do not cover the early ages of the 
evolution of stellar cl usters. In this work w e on ly present the 
result s for the auoted lBertelli et al.l (^O^) and lGirardi et all 
i200(]|) set of isochrones. 

Our second set has been computed fr om the list of 
non-rotating tracks of the Geneva group j Schalle r et alJ 
Il992l: ISchaerer et all Il993allbl: ICharbonnel et al.lll993^ ■"li> 
eluding standardmass-loss rates computed at five metallici- 
ties: Z=0.040, 0.020 (solar), 0.008, 0.004 and 0.001**. These 
tracks include all the phases of the stellar evolution from 
the zero-age main sequence to the end of the giant branch 
or the carbon ignition. The Geneva isochrones are calcu- 
lated by means of a parabolic interpolation between the 
tracks (log M, log tfc ) in the HR diagram on a variable mass 
grid as prescribed in lCervino et al.l I 200ll) . This prescription 
is also used for the isochrone computations in Starburst99. 
The isochrones were generated at specific ages with a vari- 
able time step: M = 1 Myr for 1-10 Myr, At = 5 Myr 
for 10-100 Myr, At = 100 Myr for 100 Myr to 1 Gyr, and 
At = 1 Gyr for 1-10 Gyr. Unfortunately this computation 
of the isochrones did not assure the FCT requirements (but 
it is partially implicit in the use of the parabolic interpo- 
lations), neither did it include the evolution of stars in the 
horizontal branch or later phases, which are quite relevant 
for th e resulting spectrum at ages older than IGyr^. We re- 
fer to lLancon et alj i200ll) and Vazquez & Leitherer (2004) 
for other details about the completeness of the evolutionary 
phases and their relevance in the isochrones. 

Since we are working with SSPs, we do not need to 
treat chemical evolution self-consistently in the evolutionary 

* Note that the default in Starburst99 are the " enhanced" , not 
the "standard" mass-loss rate. 

^ The Geneva group also produces evolutionary tracks from the 
begining of t he horizontal branch to t he first thermal pulse in the 
AGB phase iCharbonnel et aljIToQ^ . However, these evolution- 
ary phases have not been included due to the variations of the 
number of equivalent evolutionary points depending on the initial 
stellar mass. 



tracks. Therefore each star generation has the same metal- 
licity during the evolution. On the other hand, the chemical 
evolution in the atmosphere of the individual stars has not 
been considered in the isochrone-atmosphere model assigna- 
tion. 

Because the isochrone metallicities do not coincide to- 
tally with the stellar atmospheres metallicities, we compute 
the models combining the isochrones and atmospheres as 
follow: isochrones at Z= 0.040 and 0.030 with the twice so- 
lar atmospheres; isochrones at Z= 0.020 and 0.019 with the 
solar atmospheres; isochrones at Z= 0.008 and 0.004 with 
the half solar atmospheres; and isochrones at Z= 0.001 with 
the 1/10 solar atmospheres. 



4 RESULTS 

We are presenting the model predictions for the spectral 
evolution of single stellar populations. First, we describe the 
spectral energy distribution in the 3000-7000 A range. Then, 
we present the spectral evolution of the hydrogen and helium 
lines (section 4.2) and of the metallic lines (section 4.3). The 
effect of metallicity on the spectra is illustrated, and the 
models are discussed for two different isochrone sets (Geneva 
and Padova). 

4.1 Spectral evolution 

Figure 2 shows the spectral evolution of a single stellar pop- 
ulation with solar metallicity at seven different ages of 4 
Myr, 10 Myr, 50 Myr, 200 Myr, 1 Gyr, 3 Gyr, and 10 Gyr. 
The first effect to notice is the decrease with time of the 
stellar population luminosity, due to the most massive stars 
dying off. The second effect is the change in the shape of 
the continuum with time, up to a few Gyr; this is a conse- 
quence of the large variation of the effective temperature of 
massive and intermediate mass stars during their evolution. 
However, after 3 Gyr, the shape of the optical continuum 
is almost constant because low mass stars evolve within a 
small range of effective temperature. This results from the 
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strong concentration of the isochrones at ages older than 3- 
4 Gyr (see Figure 1). Additionally, Figure 2 illustrates the 
strong time evolution of the Balmer and the 4000 A breaks. 
The evolution of the strength of the Balmer lines is also 
striking. They have a maximum around 400 Myr, when A 
stars dominate the optical continuum. 

Figure 2 also shows that Geneva and Padova models 
have similar overall spectral characteristics; however, several 
detailed differences are clearly visible in this figure. First, the 
spectra obtained with the Padova models are always more 
luminous than the corresponding Geneva spectra. Second, 
the continuum shapes diverge in several stages of the evo- 
lution. For example, at older ages (^2 Gyr), the Padova 
spectra are systematically redder than Geneva's due to the 
inclusion of more evolved phases after the tip of the red 
giant branch, and the absence of low-mass stars^". These 
phases and low-mass stellar tracks are not included in the 
Geneva isochrones. However, at younger ages (sC 100 Myr), 
there is no systematic change of the shape of the contin- 
uum, with the Padova spectra bluer or redder than Geneva 
depending on the given age. These differences result directly 
from the isochrones (Figure 3). For example, at 50 Myr, the 
main sequence turnoff in the Padova isochrones is bluer than 
in Geneva, and the coolest stars in the Geneva tracks have 
lower TbS than in the Padova models; these differences make 
the Padova spectrum bluer than the Geneva spectrum at 50 
Myr. A significant difference between both sets of models is 
at RSG ages, i.e., 10 Myr at solar metallicity. At this age, 
the Geneva spectrum is steeper and shows stronger metallic 
bands at red wavelengths (A ^ 5500 A) than the Padova 
spectrum. These differences are produced by a larger contri- 
bution of cool and low gravity stars (Teff ~4000 K, and log 
g ~ 0) to the total luminosity in the Geneva with respect to 
the Padova models (Figure 3). 

The metallicity effect on the spectral evolution is shown 
in Figure 4. We plot the spectra of a 1 Gyr stellar popula- 
tion at four different metallicities. The effect is appreciated 
in the slope of the continuum and in the strength of the 
lines. At this age, lower metallicity spectra are bluer, have 
stronger Balmer lines, and weaker metallic lines (Figure 4a). 
These changes are produced by the shift of the isochrones 
to higher temperatures with decreasing metallicity, which is 
a consequence of stars being brighter and hotter at lower 
metallicity. The effect of metallicity is also evident in spec- 
tra at older ages (>few Gyr), where the continuum shape is 
similar in all the models, but there is a strengthening of the 
metallic lines from the most metal-poor to the most metal- 
rich stellar populations (Figure 4b) . 

4.2 Helium lines 

The high spectral resolution of our stellar library allows the 
prediction of the strength of many weak helium absorption 
lines that form in the optical continuum of young stellar 
populations. In Figure 5, we have labeled some of the most 
relevant Hei and Heii lines in the spectrum of a 7 Myr stellar 

-"^ Note that the Geneva tracks are truncated at 0.8 M0 while the 
Padova isochrones include stars with mass 0.2 M0. According 
with Vazquez & Leitherer (2004), the differences between Padova 
and Geneva models are mainly due to this truncation. 
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Figure 5. Spectrum predicted for a 7 Myr instantaneous burst at 
Z=0.001. Hydrogen and helium lines in the 3650-5000 A interval 
are indicated. 

population at Z=0.001. They are: Hel A3819, A4026, A4387, 
A4471, A4922, and A5876; and Hen A4200, A4541, A4686 
and A5412. Hen lines are detected only during the first few 
years of the evolution; their equivalent widths decreasing 
with time from values of about 0.5 A. Hel lines are more 
prominent and are detected during the first 100 Myr of evo- 
lution, when the optical continuum is dominated by B stars. 
Some of these Hel lines have been calculated previously by 
us (GLH99). 

We measured the equivalent widths of three of these 
lines (Hel A4026, A4471 and A5876) in windows of 9-14 A, 
integrating the fiux from the pseudo-continuum obtained by 
a first-order polynomial fit to the continuum windows de- 
fined in Table 1. Figure 6 shows the equivalent widths of 
these lines for a stellar population at solar metallicity. Hel 
A4026 is one of the stronger lines, while Hel A5876 is much 
weaker. The Hel lines present a maximum value between 20 
Myr and 50 Myr when the stellar population is dominated 
by B stars. Hel lines do not form after 100 Myr, the lifetime 
of B stars. The equivalent width has also a maximum be- 
tween 7-15 Myr due to the presence of post-main sequence 
stars with T^ft ^ 8000 K (see Figure 7 in GLH99). 

Note, however, that these predictions for hot stars 
are estimated using static stellar atmosphere models 
jLanz fc Hubenv^l2003^ . They do not account for wind ef- 
fects, which are important in very massive and evolved O 
stars. The lines most affected by winds are Ha and Hell 
A4686. The wind effects are less important in the Balmer 
and H el lines in the blue s pectral range, and at low metal- 
licity iWalborn et alJll995^ . 

The main effect of the metallicity on the Hel lines re- 
sults from the effect of the metallicity on the temperature 
of the turnoff, with metal-poor stellar populations having 
a higher turnoff temperature. The equivalent width of Hel 
lines is a function of Tcff and gravity, with a maximum at 
20000 K and logp = 5.0 (GL99). Thus, the variation of 
EW(Hel) with metallicity depends on the variation of the 
turnoff temperature of the stellar population with metallic- 
ity. 

4.3 Hydrogen lines 

Figure 7 shows the synthetic profiles of the Balmer lines (Ha, 
H/3, H7, H(5, He, H8, H9, HIO,..., to the Balmer limit) for 
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Figure 2. Spectral evolution for a single stellar population model at solar metallicity following the Geneva (left) evolutionary tracks 
with standard mass loss, and the Padova (right) tracks. Ages from top to bottom are 4 Myr, 10 Myr, 50 Myr, 200 Myr, 1 Gyr, 3 Gyr, 
and 10 Gyr. 



Table 1. Hydrogen and helium lines and continuum windows 



Line 


Line windows 


Continuum window (blue) 


Continuum window (red) 


H(5 


4092-4112; 4087-4117; 4070-4130 


4012-4020 


4158-4169 


H7 


4330-4350; 4325-4355; 4310-4370 


4262-4270 


4445-4453 


H/3 


4852-4872; 4847-4877; 4832-4892 


4770-4782 


4942-4954 


Ha 


6553-6573; 6548-6578 


6506-6514 


6612-6620 


Hel A4026 


4020-4031 


4012-4020 


4158-4169 


Hel A4471 


4464-4478 


4665-4675 


4464-4478 


Hei A5876 


5871-5880 


5835-5845 


5904-5912 



stellar populations at three different ages, 10 Myr, 100 Myr 
and 500 Myr at half solar metallicity, and for 1 Gyr at several 
metallicities (Z=0.019, 0.008 and 0.001). As it is well known, 
the profiles of the Balmer lines change significantly with age. 
The tendency is that the width of the lines increases with the 
evolution until about 500 Myr if the stellar population has 
solar metallicity. At ages younger than 1 Gyr, the profiles of 
these lines have only a small dependence on the metallicity. 
In individual stars, the strength of the hydrogen lines is 
independent of metal content. Therefore, the lines depend 
on metallicity only through the dependence of the stellar 
evolution in the integrated light of a stellar population. 

Since at low metallicity a star of a given mass is hot- 
ter and evolves more slowly in the HR diagram, the maxi- 
mum strength of the Balmer lines occurs at older ages. How- 
ever, the dependence of the indexes of the Balmer lines with 
metallicity is stronger (Figure 8) at ages older than 1 Gyr. 
So, a metal-poor old stellar population can have Balmer lines 



similar to an intermediate age stellar population at solar or 
over-solar. 

Figure 8 shows the evolution of the equivalent widths 
of the Balmer lines (Ha, H/3, H7 and H5). The equivalent 
widths of these lines are measured in windows of width 20 A, 
30 A, and 60 A, centered on the lines, except for Ha which 
is measured in only two windows. However, only the values 
corresponding to the window of 30 A are plotted in Figure 
8. These windows have been chosen to evaluate the effect of 
line broadening on the EW. A narrow window is required 
to date old stellar populations, because it includes less con- 
tribution of neighboring metallic lines; in contrast, a wide 
window that includes the entire line wings is more useful to 
date intermediate age populations. The measurements are 
obtained by integrating the flux pseudo-continuum found 
through a first-order polynomial fit to the continuum win- 
dows defined in Table 1. H/3, H7 and H(5 have a very similar 
behavior if the stellar population is younger than 1 Gyr. 
Their strengths range from 3 to 12 A. During the first few 
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Figure 3. Upper panel: Comparison of the isochrones obtained with the Padova (black line) and Geneva (red dotted line) tracks for 
solar metallicity and ages of 50 Myr, 200 Myr and 2 Gyr. Lower panel: Gomparison of the Geneva (left) and Padova (right) isochrones 
at 10 Myr weighted by the number of stars in each luminosity bin. Note that the relative contribution to the total luminosity of stars 
with 'Tf^ff= 4000 K is significantly larger in the Geneva than in the Padova models. 
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Figure 4. Spectra predicted for a single stellar population at different ages for different metallicities (as labeled), obtained with the 
Geneva (left) and Padova (right) isochrones. 



million years, the EW is almost constant at ~ 3 A. Then, the 
equivalent width increases with time until 400 Myr, when 
stars of 10000-9000 K dominate the optical continuum. The 
time evolution of Hq is smoother, changing from 3 to 6 A. 
Note that the evolution of the strength of these lines is very 
similar in the Geneva and Padova spectra if the stellar pop- 
ulation is younger than 1 Gyr. However, Padova models will 



predict younger ages than the Geneva models in old stellar 
populations. 

We have not plotted the equivalent widths of the high- 
order Balmer series. They have a similar behavior as H/3 and 
H7, with values that range from 2-3 to 10 A (GLH99). These 
lines are very useful for constraining the stellar ages in the 
nuclei of galaxies with an important contribution from nebu- 
lar emission (e.g. iGonzalez Deleado. Heckman. fc LeithereJ 
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Figure 7. Synthetic spectra from 3700 to 6700 A predicted for SSPs at ages 10 Myr (blue), 100 Myr (red), and 500 Myr (black), 
computed with the Padova isochrones at metallicity Z=0.008. Lower-right panel: H7 at 1 Gyr and at three metallicities, Z=0.019 (red), 
0.008 (black) and 0.001 (blue). 



l20nifl due to the strong decrease of the nebular emission with 
the Balmer decrement, while the EW of absorption lines is 
constant or increases with wavelength. 

4.4 The 4000 A break 

The Balmer discontinuity is another important age diagnos- 
tic in the optical continuum of a stellar population. Figure 
2 shows how this discontinuity evolves with time up to con- 
verging with the 4000 A break. While the 4000 A break is 
produced by the accumulation of metallic lines around 4000- 
3800 A, the Balmer discontinuity depends on the accumula- 
tion of the Balmer series below 4000 A. Different definitions 
can be adopted to measure this discontinuity, and the use- 



fulness of the definition depends on the stellar age that we 
are interest ed in measuring. In section 4.5, we will adopt the 
iRosj dlQQ'dl definition for measuring the Balmer discontinu- 
ity to date young stellar populations. H owever, we adopt the 
4000 A break definition introduced by iBaloeh et''al.l (Il999l) 
to study the general change of the discontinuity for inter- 
mediate and old stellar populations. Balogh et al. use 100 A 
continuum bands (3850-3950 A and 4000-4100 A) to mea- 
sure the break. The main advantage of this definition is that 
the index is less sensitive to reddening effects. Other defi- 
nit ions using broa der c ontinuum bands hav e been proposed 
bv lBruzuai jl983l) and lGorgas et alJ l|l999ri . 

Figure 9 shows the strong variation of the index for 
ages older than 100 Myr. This index displays a similar be- 
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Figure 8. 1) Lower panels: Equivalent widths of H(5 (left), H7 (central) and H/3 (right) at different metallicities, Z= 0.020 or 0.019 
(stars), Z= 0.040 (circles), Z= 0.008 (squares), and Z= 0.001 (crosses) from the Geneva tracks. 2) Upper left panel: Equivalent widths of 
H7 at solar metallicity from the Geneva (stars) and Padova (circles). 3) Central and right upper panels: Equivalent widths of the Balmer 
lines, H« (squares), H/3 (crosses), H7 (circles) and H(5 (stars), predicted from the Padova and Geneva tracks at solar metallicity. 



havior for the Padova and Geneva stellar models. However, 
the predictions differ for stellar populations older than 1 Ga; 
Geneva models, which do not contain all the relevant evo- 
lutionary phases, predict older ages, similarly to the case of 
the Balmer lines. On the other hand, for ages older than 
400 Myr, this index definition significantly depends on the 
metallicity, mainly due to the fact that the blue continuum 
window includes the Call K line which starts to be promi- 
nent at intermediate and old ages. The blue window also 
includes a CN band, which is strong in metal-rich old stellar 
populations. 



Recently, the 4000 A break has been used in combi- 
nation with the H5 absorption line to date stellar popula- 
tions of galaxies (Kauffmann ct al. 2003a bl3). In that work, 
the au thors use the H5a definition bv liWorthev fc Ottavianil 
1^993). This index uses a central (4082-4122 A) bandpass 
bracketed by two pseudo-contiimum band passes (4030-4082 
A and 4122-4170 A). The evolution of the 4000 A break and 
the H5a index are plotted in Figure 9. Padova and Geneva 
models show similar values in the D4ooo-H(5yi plane for ages 
younger than 2 Gyr, but, again, the Geneva models predict 
older ages. This plot also shows the metallicity dependence 
of these two indexes. For a stellar population younger than 
about 400 Myr, the D4ooo-H5yi plane does not depend on 
metallicity, but for older ages, the two indexes change sig- 
nificantly as a consequence of the variation of the 4000 A 
break with the metal content and the effect of the metals on 
stellar evolution. 



4.5 Spectral indexes 

One of the most important problems in interpreting ob- 
served galaxy spectra is related to our ability to disentangle 
age and metallicity effects. It is well known that the best way 
to decouple these two effects is by combining the strength 
of lines which are mainly sensitive to the main sequence 
turnoff temperature, like the Balmer lines, with metallic 
lines which are sensitive to the tem perature of the giant 
branch, which depends on metallicity |\yorthev et alJll994l : 
iBuzzoni. Mantegazza. fc Gariboldilll994^. T he Lick system 
jBurstein et al.l Il984l: iGorgas et al.l Il993l : IWorthev et"all 
ll994l:IWorthev fc Ottavian^l997^ . which includes more than 
25 spectral indexes at optical wavelengths, was designed to 
break this age-metallicity degeneracy. These indexes are cal- 
ibrated using Galactic stellar spectra observed at a resolu- 
tion of ~ 9 A, and are useful to date stellar populations older 
than 1 Gyr. The lack of a significant number of high and in- 
termediate mass stars in the stellar catalogs limits the ability 
of the system to date young and intermediate age popula- 
tions. However, owing to the high-spectral resolution of our 
new stellar library, the evolutionary models presented here 
are rather useful for investigating an optimum set of spec- 
tral indexes to disentangle the metallicity and age effects 
in a large range of a ges. This study will be presented else- 
where. Very recently I Jimenez et all (|200^ have pointed out 
that the best way to break the degeneracy is by including 
all the spectral features blue-wards of 4500 A. 

Alternatively, iRosel (Il984l . has proposed to use 

the ratio of the measure of the central line depth of two 
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Figure 9. (1) Evolution of the 4000 A break (D4000) measured in the synthetic spectra (upper-left panel) at solar (stars), and half solar 
(squares) metallicity for Geneva (black and green) and Padova (red and blue) models. (2) Strength of H(5^ as a function of the 4000 A 
break for different stellar models and different metallicities, Z= 0.020 or 0.019 (stars), Z= 0.040 (circles), Z= 0.008 (squares), Z= 0.004 
or Z= 0.001 (crosses), for ages between 10 Myr and 10 Gyr. 



neighboring spectral stellar features. This index definition 
has several advantages with respect to an equivalent width 
index. In particular, it is independent of the uncertainty as- 
sociated with the location of the pseudo-continuum. In addi- 
tion, this index is insensitive to reddening, but it depends on 
the spectral resolution (Vazdekis 2001; Martins et al 2004). 
Many of these indexes have been used to constrain the stellar 
population in elliptical galaxies (e.g. Rose 19 94j), but also in 
youn g and intermediate age populations lILeonardi fc R.osa 
|2003). Two relevant indexes are H(S/A4045 and Call. The 
former is defined as the ratio of the residual central inten- 
sity of RS to that of the Fel A4045 line. The latter is defined 
as the ratio of the Call H + He and Call K lines. These two 
indexes combined with the Balmer discontinuity, defined as 
the ratio of the average fi ux in two narrow bands at 3700 - 
3825 A and 3525-3600 A iRose. Stetson, fc TripiccolllQSTD . 
can constrain the age and metallicity of a stellar popula- 
tion. The behavior of these indexes with age and metallicity 
is different. H(5/A4045 is very sensitive to age, and it shows 
a strong evolution from few 100 Myr to 10 Gyr (Figure 10) , 
increasing with age. The minimum value occurs when the 
continuum is dominated by 10000 K stars, due to an en- 
hancement of H(5 in early A stars. Because stars of type 
earlier than AO have shallower Balmer lines, the index in- 
creases again toward younger ages. The Call index shows a 
similar behavior to the H(5/A4045 index at earlier ages, when 
the Call lines are weak and the ratio depends mainly on the 
He depth. Then, the index remains almost constant when 
the Balmer lines (age^ 1 Gyr) are weak and the core of the 
Call lines are saturated (this happens in stars cooler than 



F5). The Call index is also more dependent on metallicity 
than the H5 and Balmer discontinuity. Figure 10 also shows 
the behavior of these indexes for two metallicities and for 
the Geneva and Padova models. 



5 COMPARISON WITH PREVIOUS WORK 

5.1 Low spectral resolution models 

As pointed out in section 2, the new high spectral resolution 
stellar library presents two important changes with respect 
to previous work, related to the inclusion of line-blanketed 
models for hot (T^ff ^ 27500) and cool (T^g < 4500) stars. 
To evaluate these improvements and to check the changes in 
the shape of the continuum, we have compared the synthetic 
high resolution spectra with low-resolution models obtained 
by Starburst99 and GALAXEV. Figures 11-13 show the com- 
parison at three different ages, 1 Myr, 7 Myr and 5 Gyr. 

At 1 Myr (Figure 11), the optical continuum at A ^3750 
A follows quite well the contirmum, obtained with Star- 
burst99 using Lejeune (Lcicunc, Buscr, & Cuisincr 1997), 
or Kurucz, stellar atmospheres. In the near-UV, however, 
the flux is lower, and it is comparable to Starburst99 -I- 
UCL stellar atmospheres iSmith. Norris. fc Crowther. P.aI 
l2002h . UCL models include O and W-R atmospheres with 
line blanketing by [Pauldrach. Hoffmam i. & Lcimon (20()^ 
and iHiUier fc Milled lll998h . This difference between the 
fluxes at near-UV is clearly illustraded in Figure 13 of 
iLanz fc Hubenvl i2003h . They show that O stars with 
Teff^ 35000 K have lower emission in the 3000-4000 
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Figure 10. Strength of Rose's index (Call, H(5/Fel A4045) and of the Balmer Discontinuity (BD) as a function of age and metaUicity. 
The indexes are measured in the synthetic spectra that follow the Padova (black) and Geneva (red) models, at solar (left panel), at half 
solar (central panel) metallicities. In the right panels, the indexes plotted are from the Padova models, at Z= 0.019 (black, stars) and 
Z= 0.008 (green, squares). Note that after 2-3 Gyr, the Call index is constant and it shows a very small dependence with the metaUicity 
(see also Vazdekis 1999). 



A interval than the Kurucz models. At older ages (not 
shown), the agreement between Sed@ (Geneva) and Star- 
burst99+Lejeune is remarkable. 

Spectra corresponding to old stellar populations ob- 
tained with the Padova isochrones are also in very good 
agreement with the gal axev low-reso lution models which 
were computed using the lKuruc3 (119911) stellar atmospheres 
models (Figure 12). 

However, there is a significant disagreement between the 
high resolution models and low-resolution models at some 
specific ages when there is an important contribution of cool 
red supergiants, for example at 7-15 Myr in the Geneva, and 
at 9-14 Myr in the Padova models, both at solar metaUic- 
ity (Figure 13). At these ages, the contribution to the to- 
tal luminosity of low gravity and cool (Teff^ 4000 K) stars 
(see the HR diagram in Figure 13) is very high. While the 
low-resolution models from Starburst99^'^ also use Kurucz 
and/or Leujene stellar atmosphere models, here we are us- 
ing PHOENIX models for these cool stars. Even though 
PHOENIX and Leujene atmospheres at low temperatures 
agree quite well for main sequence stars, PHOENIX pro- 
duces less flux than Kurucz at red wavelengths in stars of 
low gravity. This is the reason why the high-resolution mod- 
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Figure 11. Comparison of the 1 Myr (Geneva) model at so- 
lar metaUicity (black, dotted line) with the low-resolution Star- 
burst99 models obtained with the stellar atmospheres of UCL 
(Smith et al. 2002) for massive and Wolf-Rayet stars (blue, con- 
tinuum line) and Lejeune (red, dashed line). 

els have a steeper slo pe than the low-resolution models. See 
iMartins et al.l (|22q3) for further discussion. 



Note, however, that this discrepancy also appears if we com- 
pare these high-resolution models with Sed@ low-resolution mod- 
els computed with the Geneva isochrones and Lejeune stellar at- 
mospheres. 



5.2 Intermediate spectral resolution models 

Next, we compare our synthetic spectra w ith th e available 
intermediate spectral resolution models bv I Vazd ekis (199^) 
and GALAXEV. These comparisons allow a check in terms of 
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Figure 13. Comparison of the high (black-points) and low (red-continuum line) spectral resolution models for a 7 Myr instantaneous 
burst at solar metallicity (left). The models are computed with the Geneva tracks using Starburst99-|-Lejeune stellar atmospheres (red- 
continuum line) and with Sed@ using our stellar library that includes PHOENIX atmospheres for Teff^ 4500 K (black-points). Geneva 
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Figure 12. Comparison of the 5 Gyr (Padova) model at solar 
metallicity with the low resolution model from GALAXEV (BC). 
The difference between the two models is shown as points around 
the line at y=0. 

the continuum shape but also in terms of the stellar line 
stre ngths. 

IVazdeki^ lll999l) computed single stellar population 
models at a resolution of ^ 1.8 A in two reduced spectral 
regions 3850-4500 A and 4800-5450 A. He uses a subsam- 
ple of ~ 500 stars from the original Ijone d il997ft empirical 
stellar library. The models are computed at several metal- 
licities and for ages older than a few 100 Myr. In Figure 14, 
we compare our and Vazdekis's models for a stellar popu- 
lation of 2 Gyr at solar metallicit y in the two wavelength 
regions covered bv lVazdekisI il999t) . The agreement between 
the continuum shape and the strength of the lines in the two 
models is exc ellent. 

Recentlv. iBruzual fc Charloli (|200^ presented 3 A res- 
olution models covering from 3200 A to 9500 A (galaxev 
code results), a spectral range similar to our mod- 
els. These predictions are based on the STELIB library 
iLe B orgne et al. 2003), which contains 249 stellar spectra 
at several metallicities ( — 2.0 ^ [Fe/H] -1-0.5), and spectral 
types from O to M9. However, the coverage in spectral type 
is not uniform at all metallicities. In particular, there is a sig- 
nificant lack of massive stars with 8000 ^ Tcff ^ 25000 and 



logp ^ 3.5. Figure 15 compares our and GALAXEV models 
at solar metallicity for two ages, 11 Gyr and 25 Myr. There 
is excellent agreement between both models at old ages in 
terms of continuum shape and line strengths. However, at 
younger ages, both models differ in the line strengths. The 
Balmer lines are stronger in our models, and the Balmer 
discontinuity is larger in our models. 



6 COMPARISON WITH LMC AND SMC 
CLUSTERS 

In this section, we compare the continuum shape and 
strength of the stellar lines of our models with observations 
of stellar clusters in the Magellanic Clouds (LMC and SMC). 
iLeonardi fc Ros3 ((2003) have obtained very good quality 
optical spectra of clusters in the LMC and SMC covering 
a large range in metallicity (from [Fe/H]=-|-0.01 to -1.97) 
and ages (from 40 Myr to 15.5 Gyr). The observations have 
an average spectral resolution of 3.2 A and cover the range 
- 3420 - 4750 A. These data have been kindly provided 
by Dr. Rose. Figure 16 compares the spectra of three of 
these clusters (NGC 1846, NGC 1831 and NGC 1818) with 
our models. Models and observations are flux normalized at 
4010 A. The models are also downgraded to the spectral res- 
olution of the observations and shown with a 0.5 flux bias 
for the sake of clarity. The age and metallicity reported for 
these clusters in the literature are labeled in the plot, as 
well as the values from our models. The agreement between 
observations and synthetic spectra is very good, considering 
an optimization of the flt was not attempted. 



7 LIMITATIONS OF THE MODELS IN THE 
RED SUPERGIANT PHASE 

In this section we discuss some of the limitations of the mod- 
els that are directly related with the limitations imposed by 
the stellar evolution models. One of the larger uncertainties 
in the stellar models is in the red supergiant (RSG) phase, 
due to the uncertainties that exist associated with the mass- 
loss rates and mixing (Yi 2003) . These limitations have been 
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Figure 14. Comparison of the 2 Gyr (Padova) model at solar metallicity with the Vazdekis (1999) models. The spectra are normalized 
at 4450 A (left) and at 5400 A (right). For the sake of clarity, the Vazdekis (1999) spectra are shown with a flux bias. The residual 
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Figure 15. GALAXEV models of a SSP of 25 Myr (bottom panel) and 11 Gyr (upper panel) at solar metallicity are compared with the 
high-resolution synthetic spectra presented here. These models have been degraded to 3 A resolution, which is the resolution of the 
GALAXEV models. The residual spectrum is the difference between the two models. 



previously noted by the poor match of the NIR colors pre- 
dicted by starburst models with the integrated hght of the 
observed stellar populations (Origlia et al. 199^. More re- 
cently, iMassev fc OlsenI ll2003fl have found that the Padova 
and Geneva models are not able to produce RSGs as cool 
and luminous as what is actually observed. The stellar tracks 
of 20 M0 produce RSGs of 6000-5000 K, but no RSGs 



stars in the LMC and SMC with T^// ^ 4500 K are ob- 
served. 

Even though at blue wavelengths RSGs may contribute 
less to the integrated hght of a stellar population than in 
the red range, we have found a very important effect at 
some especific ages, e.g. 5 Myr (Geneva models) and 8 Myr 
(Padova) at solar metallicity. There is a strong broad feature 
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Figure 16. Intermediate resolution spectra of three LMC clusters, NGC 1818 (bottom), NGC 1831 (center) and NGC 1846 (top), are 
compared with the models (smoothed to the resolution of the observations). The metallicity and age reported in the literature and in the 
models are indicated. The spectra are flux normalized at 4010 A. LMC spectra are shown with a 0.5 flux bias for the sake of clarity. The 
difference between the data and the model fitted is shown as a residual line around y=0. An optimization of the fit was not attempted. 



around the CallK line (at 3933 A) that is produced when 
RSGs of Te//= 6000-5000 K contribute significantly to the 
continuum (Figure 17). This peculiar feature is another clear 
indication that evolutionary tracks produce RSGs that are 
much hotter than those observed. 

At older ages (^10 Myr), high- metallicity stellar mod- 
els are able to produce RSG stars with Te// ^5 4000 K. 
However, low-metallicity stellar models also fail to repro- 
duce th e NIR colors observed in blu e compact dwarf galax- 
ies (e.g. IVazauez fc Leithereil l20o3) and stellar clu sters in 
LMC and SMC that ar e in the RSG phase (e.g. iMavval 
119971: lOrigha et aLlllflflSTl . These results suggest clearly that 
isochrones are not adequate at low metallicity and at ages 
between 10-30 Myr. For this reason, we have computed a set 
of evolutionary models that follow the stellar tracks at solar 
metallicity but where the stellar atmospheres have 1/10 and 
1/2 solar metallicity. Even though this exercise can not solve 
the problem inherent in the stellar tracks, it can provide a 
set of models that predicts some stellar characteristics that 
are more similar to those observed in stellar populations 
at low metallicity in the RSG phase. Figure 18 compares 
the spectra obtained with the Padova models that follow 
the isochrones at solar metallicity with the atmospheres at 
solar, half and 1/10 solar metallicity, and models obtained 
with the Z=0.004 stellar tracks. The main effect is seen in 



the slope of the spectra, almost identical in all the mod- 
els that follow the solar tracks, but flatter, therefore redder, 
than the spectrum that follows the tracks at Z=0.004. The 
variation of the equivalent widths of the metallic stellar lines, 
such as CallK and G band, is small at the RSG phase, and 
it is only important in stellar populations older than 1 Gyr 
(Figure 18). 



8 SUMMARY AND CONCLUSIONS 

We have computed evolutionary synthesis models that pre- 
dict high-resolution spectra of a stellar population in the 
3000-7000 A range. These models predict the indexes of 
metallic lines and the absorption-line profiles of the hydro- 
gen Balmer series and helium lines as a function of age (1 
Myr to 17 Gyr) and metallicity (Z= 0.040, 0.020, 0.008, 
0.004 and 0.001) for an instantaneous burst. The models 
are generated using the isochrone synthesis method. These 
isochrones were built based on the Geneva and Padova 
evolutionary stellar tracks. A new stellar library based on 
stellar atmospheres models is used to compute the near- 
UV and optical continuum. The library, composed of 1650 
spectra, incorporates the lastest improvements in stellar 
atmospheres: non-LTE line-blanketed models for hot stars 
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Figure 17. Spectrum (left) and isochrone (right) of a 5 Myr stellar population obtained with the Geneva isochrones at solar metallicity. 
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Figure 18. Upper panel: spectra obtained with the Padova isochrones at several metallicities: isochrone at Z= 0.020 and the stellar 
atmospheres at solar (black), half solar (green) and 1/10 solar (blue) metallicity; isochrone at Z= 0.004 and the stellar atmospheres with 
half solar metallicity (magenta). These spectra are for SSP of 14 Myr. Lower panel; equivalent widths of the CallK (3933 A) line and 
the G band (4300 A). They are measured using the line and continuum windows defined by Bica &; AUoin (1986). The colors represent 
models obtained with the isochrone and atmospheres with the metallicities indicated above. 



(Te// ^ 27500 K) and PHOENIX LTE line-blanketed mod- 
els for cool stars (3000 sST^// ^ 4500 K). Together with Ku- 
rucz models (for stars with 4750 sST^// ^ 27000 K), these 
atmospheres are used to generate high-spectral resolution 
spectra with a sampling of 0.3 A. The full set of evolutionary 
models and the stellar library spectra are available for re- 
trieval at the websites http : //www . iaa . csic. es/ ~rosa and 
|http : //www ■ iaaT csic . es/ ~mcs/Sed@ or on request from 
the authors. 

The results of the synthetic models with the Padova and 
Geneva isochrones are compared. The models obtain opti- 
mal results, or at least results as good as those of GAL AXE v 



or IVazdeki'3 il999l) for old ages. At such ages, the Padova 
isochrones work better but only in the sense that they are 
more complete than Geneva in evolutionary phases. They 
are better than previous work for young ages due to the 
larger atmospheres coverage. Geneva and Padova produce 
quite similar results. In general, age determinations are de- 
pendent on the isochrone assumptions. These models have 
limitations mainly associated with the limitations and un- 
certainties in stellar evolution. These limitations can be rel- 
evant in the RSG phase, and in particular at low-metallicity. 

Although the coverage is optimal, it is not cor- 
rect to use any model for clusters with total masses 
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below ten ti mes the Lowest Lu mi nosity Limit es- 
tablished by [Cervifio fc Luridianal i2004l): s ee als o 
ICeryino fc Valls-GabaudI 1)200,-^) : ICeryino et al.l l)2003l) . 
This limit corresponds to the situation where the integrated 
light of the cluster equals to the luminosity of a single 
star in the model, and corresponds typically to cluster 
masses around 10^ Mp), but yarie s with wayelength and 
age (see ICeryino fc Luridianal |2004| . for more details). The 
monochromatic Lowest Luminosity Limit and the intrinsic 
dispersion in the monoc hromatic flux due to incomplete 
sampling (A/" or N^ff, see lCeryifio et all (|20o3) for details) 
are included in the output of Sed® and in the tables in the 
web seryer. 
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